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ABSTRACT 


The objective of this project was to study the effect 
eaau Small artificial cavities have on natural convection 
imeon a2 HNOrizontal surface. Tests were run with water and 
Peeon il3. Data for heat flux as a function of bulk temper- 
mere difference were carefully obtained. These data yielded 
Nusselt number as a function of Grashof number or Rayleigh 
number. All of these were then compared with the data obtained 
by O'Connor [Ref. 3]. Experimental results are presented 
mometune heat transfer from horizontal circular disks, with 
mma Without artificial cavities. 

Whe data from a mirror finished surface was compared to 
that of previous investigations and the agreement was good. 
Me mwarcityeial cavities were found to affect the natural 
Poenuvection heat transfer from a horizontal surface: the 
heat transfer coefficients in both cases (water and Freon 113) 


were decreased. 





ir 


EVE, 


V. 


TABLE OF CONTENTS 


PR eneclLON==—{—{————{— =< — = 
A. BACKGROUND------------------------------ 
B. THESIS OBJECTIVES----------------------- 
EXPERIMENTAL APPARATUS—--------------------- 
A. DESIGN CONSIDERATIONS ------------------- 
B. TEST SECTION - HEATER ASSEMBLY---------- 
C, INSTRUMENTATION------------------------- 

i Temperature Measurement-—------------ 

2. Power to Heater--------------------- 
EXPERIMENTAL PROCEDURE-—--—---—----—--------.--- 
A. TEST SURFACE PREPARATION---------------- 
a Caco si ——<———<—— = 
2. Cylindrical Cavities---------------- 
TEST SURFACE INSTALLATION--------------- 
APPARATUS ASSEMBLY---------------------- 


FLUID PREPARATION----------------------- 


ef Sree Ge 


TEST PROCEDURES------------------------- 
1. Temperature------------------------- 


2. Determination of Nusselt, Grashof 
and Rayleigh Numbers---------------- 


RESULTS AND DISCUSSION---------------------- 


CONCLUSIONS AND RECOMMENDATIONS ——--—-—-—-——-—-——— 


APPENDIX A - THERMOCOUPLE CALIBRATION PROCEDURE--- 


APPENDIX B - SAMPLE CALCULATIONS ------------------ 


Noe ie er cnr=nNORRTAINTY ANALYSIS———-—-—--—---—------ 


12 
12 
18 
iy) 
at, 
29 
30 
30 
al 
33 
33 
33 
33 
34 
34 
36 
36 


36 


37 


43 
AQ 
58 
60 


68 





BIBLIOGRAPHY 


— os Se eee ee ee ee eo ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee ee es es ee ee ee ee es es ee ee ee ee 


INITIAL DISTRIBUTION LIST------------------------- 





ea 


III, 


IV. 


VI, 


VII. 


LIST OF TABLES 
Summary of Results for Water 
(Seven-Hole Test Section )------------------- 


Summary of Results for Freon 
(Seven-Hole Test Section)------------------- 


Summary of Results for Freon 
(Blank Test Section)------------------------ 


Uncertainties in the Experiment for Water--- 


Uncertainties in the Experiment for Freon 
(Seven-Hole Test Section )------------------- 


Uncertainties in the Experiment for Freon 
(Blank Test Section)------------------------ 


Uncertainty of Variables-------------------- 


40 


Al 


42 


46 


47 


48 


68 





et. 


He 


eS. 


14, 


iS. 


nO. 


Nie 


ho. 


LIST OF FIGURES 
Crossceoeetronal View of Tank, Test Section and 
Heater Assembly-------------------------------+ 


Photograph of Complete Apparatus Including 
ae ee 


Photograph of the Two Test Sections—-------- ae 
Blank Test Section--------------~~------~------~-- = 
Seven-Hole Test Section-----------~-----~-----~- 
Test Section and Heater Assembly----------~----- 
Layout of Seven-Hole Drilled Cavities---------- 
Simple Circuit for Measuring Power to Heater--- 


Log-Log Plot of Nusselt Number vs Grashof 
Number for Water With Seven—-Hoie Test Section-- 


Log-Log Plot of Nusselt Number vs Rayleigh 
Number for Water With Seven-Hole Test Section-- 


Log-Log Plot of Nusselt Number vs Grashof 
Number for Freon With Seven-Hole Test Section-- 


Log-Log Plot of Nusselt Number vs Rayleigh 
Number for Freon With Seven-Hole Test Section-- 


Log-Log Plot of Nusselt Number vs Grashof 
Number for Freon With Blank Test Section------- 


Log-Log Plot of Nusselt Number vs Rayleigh 
Number for Freon With Blank Test Section------- 


Log-Log Plot of Nusselt Number vs Grashof 
Nui pee JPee > GUI Dea OS 


Log-Log Plot of Nusselt Number vs Rayleigh 
Number for all Data—--—-——-—--—----—--------------- 


Pepemminavlolmeol ourlace Lempecrature (Water )-—- 


Determination of Surface Temperature (Freon)--- 


21 


22 
23 
24 
29 
26 
2d 


oy 


20 


SE 


O24 


28) 


o4 


20 


56 


o”¢ 
oul 


64 





SYMBOL 


D 


g 


SS 


6x 


NOMENCLATURE 


DESCRIPTION 


Heater Surface Diameter 
Aceeleration Of Gravity 


Convection Heat Transfer 
Coefficient 


Height of Cylindrical Tank 


Thermal Conductivity of 304 
Stainless Steel 


Thermal Conductivity of Test- 
Fluid Distilled Water or 
Freon 113 


Radius of Cylindrical Tank 
Bulk Temperature of Fluid 
Test Section Surface Temperature 


Temperature Difference Between Test 
Surface and Fluid (Te - Ty? 


Temperature Measured in Stainless 
Steel Test Section at a Location 
0.227 inch from Test Surface 


Temperature Measured in Stainless 
Steel Test Section 0.750 inch from 
Test Surface 


Temperature Difference in Test 
Section (T, - De 


Film Temperature Ts " ap 


Distance in Stainless Steel 
section Between Thermocouple 
Now lt and Thermocouple No. 0 


UNITS 


Inches or ft 


Ft/ a0? 


BTU /Hr/Ft? /OF 


Inches or ft 


BTU /Hr /Ft/°F 


BTU /Hr /Ft/OF 


Pnehes or 1t 


OF 


Op 


ry 


Inches or {ft 





6 el t Uaieertaimty in Variables Denoted Units of 

P by Subscript, for Example: W T variables denoted 

ie SS by subscript 
Ss 

a Pita inerna! Ditiusivic y Ft7/sec 
8 Fluid Coefficient of Thermal 

Expansion ily as 
Vv Fluid Kinematic Viscosity Ft7/sec 
DIMENSIONLESS NUMBERS 
Nu Nusselt Number Based on Test 

section Diameter hD/, 
Pr Fluid Praudtl Number le 
cr Grashof Number Based on Test Bp 6T 

Surface Diameter aia a 

v 
Gr Grashof Number Based on 3 
e Cylindrical Tank Radius BER ST 
Vv 

Ra Rayleigh Number Gr? 
p Power into Heater Watt 


10 





ACKNOWLEDGMENT 


theswaurenor wishes to express his appreciation to Dr. 
lmeeeekellenhner, Professor of Mechanical Engineering, Naval 
Postgraduate School, for his assistance, encouragement and 
Pmemunselfiish giving of his time. Without his perserverance 
this work could not have been completed. 

He would like to thank Mr. Tom Christian, Mr. George 
Bixler, Mr. Jack McKay, Mr. Roy Edwards and all of the 
Mechanical Engineering Department for their advice and work 
on design of the experimental apparatus. 

The author also expresses his gratitude to Professor 
Paul eeseos Professor Thomas Copper, Professor Dong Nguyen 
mon their advice and assistance. 

A special word of thanks and sincere appreciation to 
his fiancee, Kim-Thoa, for a long time and encouragement 


waiting. Thank God for her presence with him. 


ide 





oe NCRODUCT TON 


A. BACKGROUND 

pome recent experimental data with liquid nitrogen 
Showed that drilled cavities in a horizontal disk can improve 
tae heat transfer coeffieient during free convection [Ref. 1]. 
These cavities were approximately 0.015 inches in diameter 
and 0.025 inches deep, and were placed in a copper mirror 
finished, surface. Tests were run with none, one, seven 
and thirteen cavities and the heat transfer coefficient 
continued to improve as the number of cavities were increased. 

Free convection heat transfer occurs when a body is 
placed ina fluid at a higher or lower temperature than that 
of the body. The difference in temperature implies that heat 
flows between the fluid and the body and also causes a change 
in the density of the fluid layers in the vicinity of the 
mwirrace, On account of thermal expansion, the fluid at the 
bottom will be lighter than the fluid at the top. Because of 
this there will be a tendency on the part of the fluid to 
redistribute itself, however, this natural tendency on the 
part of the fluid will be inhibited by its own viscosity. 
The difference in density leads to downward flow of the 
heavier fluid and upward flow of the lighter. In both free 
mimnoreed CONVection, the currents transfer internal energy 
stored in the fluid, however, the intensity of the mixing 


motion is generally greater in forced convection so the heat 
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@ransicr coeificients are greater than in free convection. 
Bae heat transier coefficients in free convection are 
Felatively low but good enough for cooling many devices in 
the electrical engineering field: transformers, rectifiers, 
transmission line, electrically heated wires as in incan- 
descent lamps and the heating elements of an electric 
furnace. Free convection heat transfer plays an important 
role in steam radiators, wall of a building, especially on 
the stationary human body in a quiescent atmosphere. A 
knowledge of free convection heat transfer coefficients is 
necessary to understand and determine the heat load on air- 
conditioning or refregeration equipments. Free convection 
causes the heat losses from pipes or other heated fluids, 

te, cool ne surfaces of bodies in which heat is generated by 
fission [Ref. 2] in nuclear power applications. 

This experiment examines the effects of surface cavities 
on heat transfer from a horizontal surface with water and 
Freon 113. 

A geometrical configuration should be a first step in 
studying the natural convection in enclosures. The geometry 
and power input situations vary throughout the literature 
of studying the natural convection in enclosures. To start 
Mmpeonmexperiment, a circular test section in a cylindrical 
container, with the same size for radius of the tank and 
height of the tank was used. Heat transfer data could be 
compared with existing results, especially with O'Connor's 


fmef. 3]. 
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Hitmen weno noter across plane vertical and horizontal 
layers of gas, as well as across annular spaces, is most 
Mmmecresting from the point of view of applications. In this 
meme OL Study, we find the problem of insulation by air 
layers, as well as the heat transfer through porous materials. 
The boundary layer thickness are no longer small in compar- 
feon, with the dimensions of the air Space so that the develop- 
ment of the flow is often strongly influenced by the shape 
of the boundaries. An unstable stratification is formed if 
the heat flow is directed upward and the lower density fluid 
is located below the cooler fluid (density is higher). This 
iieeaole stratification must breakdown at a certain value of 
Rayleigh number (Gr-:-Pr) above which a convective motion must 
| be Peneraced [Ref. 4], Jeffreys [Ref. 5] and Low [Ref. 6] 
earculated this limiting value of Rayleigh number and the 
value was 1700 in case of the air layer is bounded on both 
eaaes by solid walls and a value of 1108 was found with a 
free upper surface. Schmidt and Saunders [Ref. 7] carried 
Sitean Optical investigation in a trough filled with water 
and have confirmed the first of these values very satisfac- 
torily. When the Rayleigh number is above the critical 
value, the convective motion in the fluid takes the form of 
quite regular hexagonals cells (also called Benard cells), 
the fluid moves upward in the center and returns near the 
sides. In case of a horizontal layer with circular side 
wees a CcOoneentric circular pattern of flow was obtained by 


SOberman [Ref. 8]. One method to describe the convective 
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heat transfer across a fluid layer is by an effective 
thermal conductivity Ke, which for a plane layer is defined 


by the equation [{Ref. 9]. 


K 


with Ke = thermal conductivity of a solid body which would 
admit the same heat flux as that occurring as a result of 
convection through the fluid layer of thickness 6 under the 


action of a temperature difference CT, ~ To.) (T, and T, are 


2 
the temperatures of the boundaries of the layer). Let 
K = the true thermal conductivity of the fluid then the ratio 
Ke/K describes the ralative increase in the heat transfer 
because of convection. Ke/K also represent the Nusselt 

number defined with respect to the length 6. The heat transfer 
across horizontal layers of air was measured by Mull and 

Reiher (Ref. 10] and the same measurements performed by 

Schmidt and Saunders [Ref. 7] with water (Pr = 7.1) ata 

Ra = 45000. The transition was quite sharp and could be 

easily determined. The curves for air and water do not 
coincide in a plot of Ke/K versus the Rayleigh number thus 

the true creeping motion can not be assumed, consequently Ra 
does not represent an appropriate dimensionless parameter. 

itera, correlation of Nu = C (Gr.Pr”)” US" ENO TECE oe eloakS 

Gase one obtains an approximate correlation of the data for 
mimeanad Water when m 1s put equal to 1.65. In the turbulent 
regime, the following equation is obtained 


1/3 
ee enomors)(Grepr-’ ~~) 


1S: 





eee ONOOO In the Jaminar region the 


wmalid for GriPr 
Seapuer Of the data is too great to permit a common correla- 
Guen. in the case of air the onset of turbulence appears 
less well defined. 

Torrance, et al [Ref. 11] designed an experimental 
model to study laminar flows induced by a simulated fire. 
They used air in a cylindrical enclosure with a heated disc 
in the center. Torrance, et al ([Ref. eo Te compared the 
experimental results with values determined from a numerical 
study. The Nusselt number based on a dise diameter versus 
the Grashof number was plotted. Excellent agreement between 


experimental and theoritical results were obtained. [In this 


experimental work, turbulence occurred when the Grashof 


= 


number (based on enclosure height), greater than 1.2 x Ton 


Lee and Hellman [Ref. 13] have a discussion on fire research 
in their papers. 

ieewreat transier from a horizontal circular disc, with 
Peer rOor finish in a cylindrical enclosure, was recorded by 


Duncan [Ref. 14] who obtained the following relation: 


i 2 ©. Sigs eae 


and a comparison of his results with the results obtained 
mumlorrance, et al [Ref. 12] theoritical and experimental 
was presented. 

In 1968, P. J. Marto, J. A. Moulson, M. D. Maynard [Ref. 1] 
fasceussecanthe mucleate pool boiling of Nitrogen with different 
Peuivecemecondattons. Pool boiling heat transfer of liquid 


nitrogen from circular, 1-in diameter horizontal disc was 
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Peucred,.  1ne Surfaces were: copper and nickel mirror 
finishes, also copper surfaces roughened, grease-coated and 
Teflon-coated. Drilled cavities were cylindrical holes of 
meameter of 0.0043 and 0.015 inch and spark cut conical 
hole of 0.022 inch diameter. Surface roughness and material 
alter the nucleate-boiling curve and the Teflon-coated has 
tomy littic effect while a grease coating Significantly 
decreases the nucleate-boiling heat transfer coefficient. 
feos Experiment, the effect of artificial cavities on 
both natural-convection and nucleate boiling was determined. 
Many data points were taken in the natural_convection region. 
In 1968, Torrance and Rockett [Ref. 12] made an 
analytical study of the natural convection induced in an 
enclosure by a small list spot centrally located on the floor. 
The enclosure was a circular cylinder, vertically oriented 
with a ratio of 1.0 for height over radius. A Prandtl number 
ore 0.7 (air) was fesuned and the Grashof number was based on 
the cylindrical height and list spot temperature. The 
Eauations Of Livid flow in axisymmetric cylindrical coordinates 
were solved by numerical method with a computationally stable, 
explicit method. The computation were performed with starting 
Soudtvions from quiescent through the initial transient to 
the fully developed flow. Solutions were obtained for Grashof 
number from 4 x ile to 4 x toe The flow patterns and 
temperature fields were illustrated with graphs displaying 
sets of streamlines and isotherms. With an extension of 
Grashof number to 4 X open a periodic vortex shielding 


developed, suggestive of the onset of instability. 


ey 











The cited references and diScusSion have described what 


ftiem been cone in the area of natural convection in enclosures. 


Be THESIS OBJECTIVES 


The objectives of this work were; 


ie 


To investigate the effects of drilled cavities on 
Epewnabural Cconvectin heat transfer in water and 
Freon 113. 

To obtain natural convection heat transfer data from 
a mirror-like surface of circular geometry. The 
meootlevceWwOULC tien be TGomnpared to that of Torrance, 
et al [Ref. 12] and also to the correlations of other 
investigators, especially with the results of 


O'Connor [Ref. 3]. 


This work would serve as a basis for future studies of 


natural convection heat transfer with different sizes and 


numbers of cavities and with different fluids. 
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II. EXPERIMENTAL APPARATUS 


A. DESIGN CONSIDERATIONS 

The test surface consisted of a stainless steel cylinder 
O895 inches in diameter installed in the center of a 
cylindrical tank 9.00 inches in diameter and 4.50 inches 
high so that it was flush with the tank floor. There are 
two test surfaces: one blank and one with seven drilled 
Gavities. A controllable heat source (using resistance 
heaters) was housed in a copper block, and attached to this 
cylinder. A calibrated resistor in series with the heater 
was used to determine accurately the input power to the 
heater. Eight thermocouples, wired in series were used to 
measure the average bulk temperature of the fluid and one 
removable thermocouple was used to provide a check on this 
temperature. To record the temperature distribution of the 
stainless steel cylinder, four stainless steel sheathed 
thermocouples were equally spaced ane the axis of the 
stainless steel test section. The temperatures at these 
locations were obtained using a Dana digital voltmeter. 
The test surface temperature was obtained by linearly 
extrapolating the two nearest thermocouples to the test 
surface (see Fig. 17). 

Figure 1 is a cross sectional view of the tank, test 
section and heater assembly. Figure 2 is an overall view 


of the complete apparatus (including instrumentation). 
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Figure 3 shows a photograph of the two test section cylinders. 
The material for constructing of the tank, test section and 
heater arrangement were available from the apparatus used by 
m Connor [Ref. 3]. 

The following criteria were used to establish the 

apparatus design: 

(1) The enclosure radius was to be equal to the enclosure 
height. This ratio was convenient to work with and 
the results obtained could be directly compared to 
the experimental and theoritical results obtained 
by Torrance, et al [Ref. 12], Duncan [Ref. 14] and 
Omeonnor f[Ref. 3]. 

(2) The radius of the test surface was to be equal to 
Ce Srcnem Of Lhe herenteen the enclosure, this is 
also the value used in the work of Torrance, et al 
fiven, wae 1). 

foe tewratioc Of the depth cavity to diameter was 
approximately two, this corresponds to the value used 
byevioulson {Rer. 15]. 

Insulation was used to minimize the radial heat loss 

due to conduction as much as possible. 

The test section designs used by Mac Kenzie [Ref. 16] 
and Duncan [Ref. 14] were examined. They both used a solid 
cylinder to get one dimentional heat flow. In Duncan's 
work, he used a stainless steel cylinder with a flanged top 
forming the test surface. Duncan [Ref. 14] discussed his 
results and showed the effects of heat transfer due to this 
design. Mac Kenzie [Ref. 16] used a similar cylinder-heater 
assembly in his work. To get one dimentional heat flow in 
cylinder we had to minimize the radial conduction losses. 


The cylinder would not be silver soldered to the heater 


section as O'Connor (Ref. 3] did, but a special design as 
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Shown in (lig. 6) for a changeable cylinder test section was 
constructed. This test section-heater assembly would be 
rigidly supported at the base of the heater but also it could 
be raised or lowered by screws (see Fig. 6) to keep the test 
Surtiace at the same level with the bottom surface of the tank. 
Because the heat flow through the cylinder would be the basis 
for the calculation of the convection heat-transfer coeffi- 
cient (h), the thermal conductivity of the material used had 
to be accurately known. For this reason and because of its 
ease of machinability, 304 stainless steel was chosen for the 
best Surface. 

To minimize the heat loss by conduction through the 
wemeeom Or the tank (since this bottom would be in contact 
with the test section), a material with a very low thermal 
conductivity was chosen. Phenolic was readily available and 
it was used by Mac Kenzie [Ref. 16] in a similar apparatus, 
thus it was Known to be strong enough to resist the stresses 
applied in achieving a watertight tank arrangement. The 
stainless steel cylinder was centered in the bottom of the 
tank, To minimize the conduction area, the phenolic was 
machined to a tight fit and beveled. In the phenolic bottom 
plate (Fig. 1) a 0.0625 inch Neoprene rubber "0" ring, 
lubricated with silicone high vacuum grease was: instalied. 


Plexiglass was used for the tank top and sides. 
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B. TEST SECTION-HEATER ASSEMBLY 

There were two test surfaces used in this experiment: 
one blank and one with seven cylindrical drilled cavities. 
The two-inch stainless steel cylinder (Fig. 3) was machined 
from a one-inch diameter 304 stainless steel rod to a .895 
miem diameter. Four 0.041 inch diameter holes were drilled 
to a depth of 0.375 inch, Separated axially approximately 
Oyo inch and 90° Apert £Or thermocouple installation. The 
Stainless steel cylinder was attached to a copper block, by 
means of four screws to facilitate changing test sections 
(see Fig. 6). The copper block served as a housing for the 
heaters. To heat the test surface, five Watiow Firerod 
MaeLridge were used, installed in the copper block. A 3.375 
inch section of 2.0 inch diameter stainless steel tubing with 
meeeO inch outside and 0.625 inch inside diameter ring 
welded to one end was attached to the copper block and 
heater by two 0.75 inch machine screws, provided as a rigid 
support for the test section-heater assembly. This also 
eliminated the problem of metal to metal contact with the 
test section. In order to minimize the heat loss through 
this support, several holes were cut. 

For seven-hole test section, seven cylindrical holes 
were drilled to a depth of 0.0375 inch and a diameter of 
0.015 inch. Figure 7 is the layout of these seven holes on 
Bion LeSt surface. 

Several layers of asbestos paper were placed between the 


heater base and stainless steel ring. The stainless steel 
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moiecemnachined to a tight fit to the pipe. These rings were 
bemtead to a 12,0 inch diameter phenolic plate, which served 
aoa SUpport plate for the entire system. One end of the 
Seeneess steel pipe was fitted anto the stainless steel 
moe, ithe test section-heater assembly was enclosed in a 
Weezo inch thick steel can with a diameter of 6.25 inch. 

It iS important to minimize the heat loss from the test 
section and heater assembly. 

several layers of asbestos paper were wrapped around the 
heater assembly and a portion of the test section and secured 
with insulating tape. Several layers of Flex-Min K insulation 
were wrapped around the entire assembly, coupletely insulating 
the interior of the stainless steel can. Several layers of 
asbestos paper were put into the hole at the bottom of the 
Stainless steel can, between the stainless steel ring holding 


the pipe section and the base plate (phenolic). 


C. INSTRUMENTATION 
1. Temperature Measurement 

The accuracy of measurement of the temperatures was 
aired Tomes eexpDerimcne: SO, IsAviype [ copper-—constant, 
an thermocouples sheathed in stainless steel (sheath diameter 
diameter = 0.035 inch) were selected. The temperature dis- 
tribution in the stainless steel section was determined using 
four of these thermocouples. The temperature of the surface 
was found using an extrapolating method (see Appendix II, 


Hee 17). 
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The bulk temperature of the fluid was determined using 
eight Copper Coustantan thermocouples wired in series (four 
pairs) manufactured a Dynatech type of welder and a single 
removable thermocouple to provide a check on this temperature. 
The eight thermocouples were wired in series into an insulated 
Switching box so the reading gave an average temperature. 

The sheathed thermocouple leads were securely fastened to the 
bottom phenolic plate of the assembly and the more plexible 
thermocouple extension wire was used to extend the wire into 
the switching box. A reference junction thermocouple was 
inserted in a sealed oil-filled glass tube which was immersed 
in a crushed ice-distilled water bath. This thermocouple was 
also wired into the switching box. 

The output from the thermocouples was read on a Dana 
aortal voltmeter accurate to 0.001 milivolt. 


2 Power to Heater 


PRECISION 
HEATER (R,,) RESISTOR (R) 





Figure 8. Simple Circuit for Measuring Power to Heater 
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In order to measure the power input to the heater, 
a circuit as shown in Figure 8 was constructed. Since Re 
fomatunction of temperature, we should use a more stable 
meormscbor (R) (calibrated resistor). With this simple arrange- 
meme, a voltage across heater (Vis) and across resistance 


(Vp) were taken. The current through the circuit could be 


determined: 
I,, = VR/R (3) 
Knowing R (2,013 3 Vis and Vp» power was calculated 
irom : 
p = Vily 
and from (1): 
P = V,(V,/R) (4) 


To get the input power to heater accurately for every 
run, we needed to measure the values of Vir and Vp from (4) 


only. 
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IfI, EXPERIMENTAL PROCEDURE 


A, TEST SURFACE PREPARATION 
1. Mirror Finish 

A procedure as outlined below was used to prepare 
each test surface; 

The test surface was first sanded with emery paper 
number 1 and 3. After that, it was washed with detergent 
and warm water to remove the previous abraSive also a final 
wash with alcohol. 

In order to get a mirror finish surface the test 
surface was then wet polished on nets ometiee dl wheel covered 
with velvet. It was impregnated with one micron diamond 
dust. The test surface was placed face down on the wheel and 
a low speed of rotation was used. Usually, it took from 
fourty five minutes Ge an hour to complete the work. During 
the work, sTlooell was used to wet the surface at all time. 
After completion of the polishing, the test surface was washed 
with alcohol again and was dried with a dryer. After work was 
completed, it waS wrapped in special paper to prevent con- 
mamination and touching. 

Weeeeveeendrical Cavities 

| The seven-hole test surface was drilled before the 
fitrOoreitinish. The purpose of making a mirror like surface 
to insure that the surface conditions were not responsible 


for the heat-flux of the surface. 
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on ees were drilled using Sphinx spirec pivot drill, 
me. 01S inch diameter, The cavity depth was controlled by 
using a fesler gauge on the lathe. A ratio of approximately 
2.9 was desired for cavity depth over diameter, So the 0.015 
inch holes were drilled to a depth between 0.035 and 0.040 
mach . 

The locations of seven holes on the test surface was 
Shown in Figure 7: One in the center of the surface and the 
other six equally spaced on a circle whose center was at the 
first hole. The holes were spaced 0.165 inch apart. 

All of the final measurements were taken to determine the 
location of the center of the holes for sheathed thermocouples 
relative to the test surface and the final test diameter 
was measured. An engineering flat and a scale readable up 


to 0.0001 inch were used, 


B. TEST SURFACE INSTALLATION 

Immediately before installation of the seven-hole test 
Surface, it was put in an ultrasonic cleaner for ten minutes 
then it was washed with distilled water, rinsed with acetone 


and dried, 


C. APPARATUS ASSEMBLY 

There are two main parts in the apparatus: 1) the 
test sections and heater assembly, 2) the tank. 

The test heater assembly (section III, B) was rigidly 
bolted to the phenolic base. After finishing the insulation 


around the test section and heater assembly, the can was 
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mastalled with the insulation inside. A phenolic disc used 
weethe bottom of tne tank was fitted to the test section 
carefully. The level of the test section was adjusted from 
the bottom base of the assembly by screws. Three long 
threaded rods were used to secure the upper plate to the 
phenolic base plate, 1LDO apart, These three rods also used 
to adjust the test surface in level with with phenolic 
bottom plate. A rubber "0" ring was fitted around the test 
Surface and the hole in the bottom plate to prevent leaking. 
For water, rubber gaskets were used at the bottom and top 

of the cylindrical plexiglass to make sure that the enclosure 
did not leak. To secure the top plexiglass plate to the 
enolic bottom plate of the tank, six threaded rods were 
used, B04 apart . 

To measure the average bulk temperature of the fluid, 
four small holes were drilled through the top plexiglass of 
the tank to insert four pairs of themocouples all wired in 
series, near the tank wall. One bigger hole was used to 
pour the fluid into or out of the tank, and was also used 
as a hole to insert one removable thermocouple used to provide 
a check on this temperature. It was possible to remove the 
plexiglass top plate and also the cylindrical wall of the 
tank without disturbing the test section in order to change 
ges fluid for another run. 

It should be noted that the test surface and phenolic 
bottom plate were cleaned with acetone prior to filling the 
tit. All other surfaces in contact with the fluid were 


cleaned with distilled water. 








D, FLUID PREPARATION 
The test fluid were distilled water and Freon 113. The 
fluid could be poured into the tank through a one inch 


diameter hole. <A siphon was used to fill the tank. 


E, TEST PROCEDURE 

1. Temperature 

For water, five settings of power input to the heater 

(10, 20, 30, 40, 50 watts) and for Freon 113 seven settings 
(1, 2, 3, 4, 5, 6, 7 watts) were used. Data were taken at 
these settings of power input. At the beginning of each 
test, a lowest power input was set by the method nilecueeed 
in section II, C 2, A long time was allowed to elapse until 
the steady state was reached (usually two to three days for 
each data point). The same procedure was repeated for each 
higher power input setting. When we reached the highest power 
input to the next lower one. When the starting point had 
reached steady state a complete set of data was obtained for 
this run. 

For each setting, we took data for: 1) themocouples 
readings in the stainless steel test section, 2) the bulk 
temperature of the fluid, 3) voltages across heater and 
calibrated resistor. 

Temperature readings were usually taken every hour in the 
daytime only. When the bulk temperature and the temperature 
nearest the heater were changed less than half degree per 


hour, steady state was assumed to have been reached, 











Peele tecimination Of Nusselt, Grashof and Rayleigh Numbers 


To determine the Nusselt, Grashof and Rayleigh numbers 


the fluid properties were calculated at the film temperature 


which is the average of the bulk temperature (T,) and the 
J belie ae 
test surface temperature (TQ) 1.e. Tl; = ae The thermal 


menauctivity of 304 stainless steel (K..) was calculated 
using the data from Touloukian [{Ref. 17] at the average of 
the temperature of the two nearest tempermocouples to the 
test surface i.e, de = “oh. This corresponds to the average 
temperature of the length Ax = (X4-XQ)- 

For each run, the temperature distribution in the stain- 
less steel was found to be linear. For this reason, the 


test surface temperature could be determined by using linear 


extrapolation 


T = B,x 4 By (5) 


Using the two nearest themocouples readings to the test 
eimrace, Substituting into equation (5) and note that, at 


surface, x = O, thus: 


Ey te. (6) 


T enmePaee 2 


and 


CH se ti 
B. a heat flux at surface = aw =o ee (7) 
1 Ax (xX, - Xo) 


A simple energy balance at the surface, using Fourier's 


Law of Conduction yields: 


o¢ 











6T. 
Kener aT 


ss 6x surface ~ \epiiee (8) 


The Nusselt number was obtained from the definition: 


Nu = —— (9) 
Ke 
and 
OeenACT. - T.) = K_. a (25) (10) 
S b Ss Axx 
or 
a ee 
or - Ty? A(By - T,) (11) 
z QD 
From (9), Nu = ACB, = 1, )KE (12) 
2 4 
From (10) and (7), we thus obtained: 
Bak 
Nua a aa Pe) (13) 
if 2 4 


The Grashof number was calculated using the diameter of 


the heated test surface as the characteristic length: 


g8 DY (T, - T,) 
Gr <5 (14) 
Vv 
or 
. 8 73 
Gr = “3 DY (By - Ty) (15) 


The Rayleigh number is the product of the Grashof number 


and Praudlt number: 
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Ra = Gr - Pr (16) 


A simple calculation for water is contained in Appendix 
B, Lt and for Freon in Appendix B, 2. 
If we use the heigh or radius of the tank as a character- 


istic length then the Grashof number becomes: 


or. = Gr = (88) H? (By - T 


: (17) 
Vv 


4) 


Since the test section diameter was chosen to be one fifth 
the heigh of the tank thus we have the relation for Grashort 


numbers: 


Gre = Gr, = 125 Gr (18) 
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TABLE I 


UNCERTAINTIES IN THE EXPERIMENT (For Water) 


Power input Percent Uncertainty 
Settings (Watts) 





Nusselt Number Grashof Number Rayleigh Number 


oi. 12 6 6 6 
20.12 4 4 4 
30.93 3 8 8 
40.93 $: 2 5 
00.54 3 i 2 
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TABLE II 


UNCERTAINTIES IN THE EXPERIMENT (Freon, Seven-Hole) 


OL DR SS a ee ee ee oie 


Power input Pereent Uncertainty 
Settings (Watts) 


Nusselt Number Grashof Number Rayleigh Number 
6 8 7 
FS 


aS 


alS 
EO 


NY OO 7 fF WO WD 
oO XN Ob Oo OO WD OW 


) 8 
c 2, 
74 8 8 
8 ) 
8 g) 
g 7 


2d 
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Power input 
settings (Watts) 


NY Oo OF F*® ©§O NH) 


-12 
13 
Pls} 
213 
.28 
-ol1 


. 26 


TABLE 


Nusselt Number 


Oo nN N OO N CO 
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UNCERTAINTIES IN THE EXPERIMENT (Freon, Blank) 


Percent Uncertainty 


Grashof Number 


6 


8 
6 
7 
8 
7 
+) 











Rayleigh Number 


6 


6 
6 
7 
7 
6 
7 





IV. RESULTS AND DISCUSSION 


This investigation covers a range of Grashof number (based 
on test surface diameter) from 1 x hohe Foro: x 10 hOr Water 
with seven-hole test section and for Freon 113 a range from 


x 10° ti@ae Sx 10° with seven-hole test section and from 


5) 4 ite” LOwmcom x ee with blank test section. A range of 
Prandlt number from 2.58 to 4.65 (with seven-hole test section) 
merewatver and for Freon 113, at a range from 7.65 to 8.84 
(seven-hole test section) and from 7.52 to 8.38 with blank 
test section. A summary of the results can be found from 
iaele IV, V and VI. 

Two different surfaces were tested and twenty-nine runs 
were made. The natural convection data points were plotted 
in Figure 9 through 16 and approximated by straight lines. 

In section II, C the bulk temperature of fluid was cal- 
culated. To reach steady state for each input setting, a 
minimum of twenty-four hours was needed. Figure 9 isa 
log-log plot of the Nusselt versus the Grashof number for 
water with seven-hole test surface. For this plot, a very 
good agreement between the heating and cooling phase of 
Nusselt and Grashof number (based on test surface diameter) 
was obtained. 

A brief discussion of the result of water with seven-hole 
test section as compared with the work of O'Connor with 


blank test section [Ref. 3] is needed to have some 














@enclusions on the effect(s) of drilled cavities on heat 
transfer (ref to Figures 9, 10). The log-log plot of the 
Nusselt number vs Grashof number (based on test section 
diameter) was approximated by a Straight line which has the 
same slope as that of the least square fit of O'Connor's re- 
Pees DUT With a lower intercept. The work of Mr. O'Connor 
Peemeiseussed ihn the introduction of this experiment, Also, 
@ecomparison of the result with that of Torrance [(Ref. 11] 
is illustrated in Figure 9. A Prandlt number of 0.7 (air) 
was assumed by Torrance and the Grashof number was based on 
cylinder height and list spot temperature. 

For Freon 113 with seven-hole test section, the approxi- 
mate straight line is lower than that of the blank test 
Section but the slopes are the same in both cases. A log- 
log plot of the Nusselt number versus the Grashof number for 
seven-hole test section and blank test section is shown in 
Figure 109. 

Senne of the data points were taken during the cooling 
phase (when power inputs were decreased). These points were 
denoted by << ' >> in the Table, For cooling phase at a 
power input Pecans, to compare with that of the heating 
phase at the same Grashof number the Nusslet number is a 
Mittle bit lower. For Nusselt number versus Rayleigh number 
the data for cooling phase (correspond to the same power 
input setting) has almost the same Rayleigh number with the 
corresponding point in the upward phase but the Nusselt 


number is lower (for the case of water with seven-hole test 
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section, Figure 9); for Freon 113 the corresponding point 
in the cooling phase has almost the same Ra but the Nu is 
lower. The Nusselt number vs Rayleigh number (log-log plot) 
is plotted in Figure 10 for water with seven-hole test section 
memecompare with the result of O'Connor [Ref. 3]. Figure 12 
meee plot for the case of Freon 113. 

The heat transfer correlation which results from this 


experimental data is: 


Nu = 0,38 (aos for water with seven-hole test section 
Nu = 0.48 “ney Vo for Freon with seven-hole test section 
and Nu = 0.53 GPa) oo for Freon with blank test section. 


Torrance [Ref. 12] in his numerical solution, assumed a uniform 
fluid bulk temperature and in his experimental work, used 

an average of the container wall temperature and the floor 
temperature. He obtained a good agreement between his 
experimental and numerical method. In this thesis, an average 
bulk temperature was used (automatically obtained from the 


reading of eight thermocouples inserted in the fluid). 
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V. CONCLUSIONS AND RECOMMENDATIONS 


The results lead to the following conclusions: 


Cy find teal “arilkied cavities affect the natural 
convection heat-transfer of water and Freon 113: the 
line of log-long plot of the Nusselt number vs Grashof 
number and Rayleigh number with seven-hole test section 
is lower than that of O'Connor (with blank test section) 
for water, but the slope is almost the same. 


The heat transfer line in the log-log plot of the 
Nusselt number versus Grashof number and Rayleigh 
number for Freon 113 is higher than that of water with 
the same test section. 


For water with seven-hole test section, a correlation 
which results from this experimental data is: 

Nu = o.BaCnay = for water and Nu = 0.48(Ra)°°? 
HOmeceon 11S wiht Seven-nole Test section and 


Nu = O,82Ceaye 


ne Freon 113 with blank test section 


0 


The following recommendations are made for future 


Continue the present studies with different size and 
number of cavities (say one, thirteen, and ninety six 
cavities) and compare the results with this experimental 
work. 


Using a power setting of at least eight watts for the 
case of Freon 113 to get a boiling curve and the result 
would be compared to that of Mr. Moulson [Ref. 15]. 


Using the same test sections but different fluids to 
investigate the effects of changing fluid. 


Study the effects of different cavity geometeries 
(says cylindrical, conijieal, reentrant). 


Find a method to get the bubbles out of the tank if 
the power input is greater than 40 watts in case of 
water to get a good result on heat transfer. 


Rework this experimental stady with the flow visual- 
iZation technique and compare the results with the 
results of this study. 
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APPENDIX A 
THERMOCOUPLE CALIBRATION PROCEDURE 


The calibration of the sheathed and bulk temperature 
thermocouples was a major part of this experiment since all 
the results and discussion based on the accuracy of the 
determination of the test surface and fluid temperatures, 

All the wiring and recording instruments would be used in 
the experiment were not changed after the calibration was 
completed. 

A Rosemount type of calibration system was used with a 
constant temperature oil path. The ice junction thermocouple 
(sealed in a glass tube, filled with oil) was placed ina 
Stainless steel dewar filled with ice water and circulated 
with a pump. The nine bulk temperature thermocouples and 
the four sheathed thermocouples were suspended about three 
about three inches into the oil bath. The calibration was 
conducted over a range from go°F to 324°F and a maximum error 
for the thermometer for this temperature range was 0.005°F. 

A Platinum Resistance thermometer with a commutation bridge 
was used as a Standard. 

A comparison of the calibration results with standard 
thermocouple tables: for the sheathed thermocouples the 
result was good with deviation from the table values in a 
mamee trom —-O.012 mv to +0.005 mv and for the bulk temperature 


ita range from -0.009 mv to +0.001 mv. 
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The above results led us to use the standard tables with 
an uncertainty in the sheathed thermocouples reading of 


0.5°F and 0.5°F in the bulk temperature reading. 








APPENDIX B 
SAMPLE CALCULATIONS 


I, FOR WATER 

timese Calculation are for data recorded on 20 December 
1974 for an input of 11.12 watts and with seven-hole test 
suriace. 
1, Determination of the Surface Temperature 

fOnaerermine the temperature of the test surface a 
graphical method was used. A plot of the temperature versus 
the thermocouple distance from the test surface was made for 
each run (corresponding to a proper setting of power input), 
as shown in Figure 17 using the two nearest sheathed thermo- 
couples to the test surface. Using a linear extrapolation 
to obtain the surface temperature an error within 0.5°F 
eould be obtained. 
Pee calculation of the Nusselt Number 

Wernaa the equation (13) of Section fil, EZ to calculate 
Bhi Nusselt number : 


B DK os 
iy a (II~-A) 
K-(Bo - Ty) 


The value of Ke was obtained based on an average 


S 
temperature of the cylinder over the length Ax using the two 
mearest temperatures to the test surface in the stainless 


steel cylinder. The values of Kss were obtained using the 





Figure 17 


DETERMINATION OF SURFACE TEMPERATURE 


Power Input Setting: 11.12 Watts 


Bulk Temperature: 85.60 as 


Temperature 
ey 
140 
Is) 
120 a 
/110.15 


110 Se 
Extrapolated Surface Temperature 


100 | = 
0 2a 6 .8 TO Seo 


Distance rom teste wourtace Cinclhies ) 


G1 








@ata frem Touloukian [Ref, 17], Fluid properties were 
jb dh 


Svaluated at the film temperature i.e. Ty a Hf. 


(section III, E2) using the data from the Table VI, Appendix 
feericom Kreith [Ref. 18]. For 241.12 Watts input setting, 


the following were obtained: 


ta 120.04°F 

ee 139, 54°F 
_ O 

T, = 85,00°F 


D = 0.07458 ft. 


Ax = 0,04183 ft. 
B, = 566.14 Btu/hr-ft” 
By = oS 1 


Tss = 129,79°F 

T, = 97.88°F 

Kss = 9.099 Btu/hr-ft-°F 

K, = .363 Btu/hr-ft-"F 

gB/,2 = 111.000 x 10° 1/£t*-or 
Pr = 4,65 | 


GC Orme) 099) 
and from (II-A) Nu -RGCGROMIG=s55c0) Song 


3, Calculation of the Grashof Number 
Wemne equation (15) of section 111, EZ to calculate the 


Grashof number: 


Gr = (g8/,2) D® (B, - T,) (15) 


The values gB/,2 Were eCalCllameamenan Appendix C, from 


hreith [Ref. 18]. 
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For i1,12 watts input setting, the Grashof number is; 


ee 1100x110" (o.07458) cilanis es5460) = 1.13 x 10° 


4, Calculation of the Rayleigh Number 
Using equation (16) of section III, E2, the Rayleigh 
mumpet 1S the product of the Grashof number and Prandlt 


number: 


Ra = Gr : Pr (19) 
uSing (12) the Rayleigh number is 
6 


R, = 1.13 x 10°(4.65) = 5.25 x 10 


i FOR FREON 113 

The following is a sample calculation of the Nusselt, 
Grashof an Rayleigh numbers for data recorded on 19 March 
1975 for 1.16 watt of power input setting and with seven- 
meee test surface. 
meeepecermination of the Surface Temperature 

Using the same method in I-1, the test surface temperature 
was obtained in Figure 18. 
2. Calculation of the Nusselt Number 

Weime the quation (lI1-1) of section Til, E2 to calculate 


the Nusselt number: 


(B, )(D)(K,,) 
Nu = =o 


K-(By - Tq) Chie) 


The values of i were obtained using the data from 


Touloukian [Ref. 17]. Fluid properties were evaluated at the 


ee eae 
film temperature qT; = a Csecetona rl win2) usine the 








Females 


DEEERMINATION OF SURKACE TEMPERATURE 


Power Input Setting: 1.16 Watt 


Bulk Temperacure: 76.17 OF 
Temperature 
OF 
110 
100 
90 
ee 27 
80 Extrapolated Surface Temperature 


790 





ae .4 AG oe) iO) 1a 


DistancemimomeTest Surface (inches) 
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data from Table |] and Table 13 from Thermophysical Properties 
Seeeirigerant 113 (Ref, 19], For 1.16 watt of power input 


Setting, the following were obtained: 


P= S4m or 

T, = 86,12°F 

T, = 76.17°F 

= O Onis =. 

BL = 45.18 Btu/hr-ft* 
By = 83, 27°F 


Tss = 85.17°F 


Vo) a ee yee 


f 
Kss = 8.81 Btu/hr-ft-°F 
_= 1.592 lbm/ft-hr 
K, = 0,0432 Btu/hr-ft-°F 
Cp = .229 Btu/lbm-°F 

0 = 97.47 lbm/ft” 


The values of 8 were calculated by definition from the 


thermo-dynamic properties [Ref. 20]. 





1 
meee dp, . | Lee eee V “ov 
B= - 5 ar? 0 SS Z ie” “Gp ee aes (20) 
0 Ts _ tT, 
with 
_ “SQ aaesal 
‘f 2 


For 1,16 watt of power input setting, we have: 


-4 1/o 


8B = 8.28 x 10 R uSing a temperature difference 


= 10°F in equation (20). 


ae) 

i 

fer | 
! 
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Tine Prandit numbers were calculated by definition; 


Pr = v/a = 


|= 
| 


Rome e  wGe: (21) 
Wsane (21), the Prandlt number for 1.16 watt of power 
input setting is: 


melo (220) 


0.0432 = 8,84 


Using equation (9), we have: 


_ (45.18)(0.07458)(8.81) _ on 
(EEEnOnOIe2 (83.27 = 76017) °°!” 


See Calculation of the Grashof Number 
Using equation (15) of section III, E2 to calculate the 


Grashof number 


Gr = (g8/,2) D” (By = Ty) (22) 


Meine the units stated in the table of the thermophysical 


properties of Freon 113 [Ref, 21] and D = 0.07458 ft, we have: 


Gr = (173135.52)(8)(p/u)"(B, - T (23) 


a? 
Using this equation, the Grashof number is: 


y 
(173235.52)(8.28 x 107*) (245) (83.27 - 76.17) 


Gr 





6 


II 


Gr Corot ) x 10 
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muecarculation of the Rayleigh Number 
Weine equacion (12) of section III, EZ, the Rayleigh 


number is the product of the Grashof number and Prandlt 


number: 
Ra = Gr « Pr (24) 
Thus, 


(3.81 x 10°)(8.84) = 33.68 x 10 


Ra 
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APPENDIX C 


UNCERTAINTY ANALYSIS 


The method proposed by Kline and Mc Clintock [Ref,. 21] 
to obtain the uncertainties was used in this experimental 
work, The following table is a summary of these quantities 
including assumptions and approximations to get an estimate 


of the uncertainties for each variable. 


AS lid ayy tet 


Uncertainty of Variables 


———. 


Variable | Basis for Value Wilecertaanty ecu) 
Ts dictated by scale of graphical 0.5°F 
method 
Ty thermocouple calibration 0.5°F 
To and T; thermocouple calibration 0.5°F 
K.(H,9) assumption based on table . 
values [Ref. 18] 0.001 Btu/hr--ft-"F 
Kss(304 s.s.) assumption based on table . 
values (Ref. 17] 5% Btu/hr-ft- F 
K, (Freon 113) assumption based on table - 
values [Ref. 19] 5% Btu/hr-ft- F 
gB/,2(H,9) assumption based on table a6 
values [Ref. 18] OVO yt tae 
gB/,2 assumption based on table ae 
(Freon 113) values [Ref. 19] 0-01 hj tte 
6x accuracy of measurement 0.42 x 10°" ft 
D accuracy of measurement 0.20 x 10°° rt 
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Pr(Ho0) assumption based on table 
values [Ref. 18] 5% 


u(Freon 113) assumption based on table 
values [{Ref. 19] 0.002 lbm/ft-—hr 


Cp(Freon 113) assumption based on table 5 
values [{Ref. 19] 0,001 Btu/lbm-~F 


B( Freon 113) assumption based on table _4 
values [Ref. 19] On xo One lo 


o0(Freon 113) assumption based on table 3 
values [Ref. 19] 0.01 lbm/ft 


Using the second-power equation of Kline and Mc Clintock 


Ref, 21 to calculate the uncertainties in the values obtained 


experimentally. 
Let Ung = uncertainty in Ts 
Ump = uncertainty in ie 
x a 
then, U,.,. = U See iy le 
me or Ts - Tb Ts Tb 
4 
= | (0.5)? 3 (0.5)? | = Fern (25) 
1 
and Us _ Un, . _ cu)? : (Ur. 7 2 
See 1 9 e 1 0 
4 
= [co.5)? £ (O85 | = 0.707°F (26) 


From equation (8), section III, E2, we have 


and using the second-power equation of Kline and Mc Clintock 


[Ref. 21] we have: 
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h Kss OX ol: Cay) 
SS 
mimeot the uncertainties are found from Table VII, 
By definition, the Nusselt number is: 
ho = aa 
f 
and the second power equation for uncertainties of this 
equation is: 
ee 
Z 
UNu = 4) : i, ee + [ Kz) C25) 
Nu h D K, 
U y) Uy . 2 
The values of | 2 | and =| can be neglected. 
D K 
f 
Thus, equation (21) becomes: 
UNu = vh (29) 
Nu h 
Using equation (15), section III, EZ, to determine the 
uncertainty in the calculation of the Grashof number: 
Gr =(g8/.2) D® (T_ - T,) 
Vv S b 
and i 2 5 9. 2 
nor [7/2] [2] 4-89] (30) 
Gr gB/,,2 D OT 


Hemerreon 113, using equation (17), section II, A3, to 
determine the uncertainty in calculation of the Grashof 


number: 
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Gr = (173135.52)(B)(p/v)* 67 


and 
2 Dy 2 2, 


i. U 2U 2U U 
Yor = (28) + (29) «(9 + GE] we 


All of the uncertainties are found from Table VII, 
Using equation (16), section III, E2, to determine the 


uncertainty in the calculation of the Rayleigh number: 
Ra = Gr + Pr 
ana 
vRa -| [-S*] 
Ra Gr 


For Freon 113, by definition 


NH 


2 Up. 2 
+(e] (for H,0) (32) 


S 
f 


ie 
.S 


Pr = 


~ 


and uncertainty in the calculation of the Prandlt is: 


mera) lop) + Leet) | 


All of the uncertainties are found from Table VII, 


Thus, for Freon 113, we have 


Nl 


U U U 
ery .c}. +{ 3H 


A summary of the results of this analysis is in Table I 
for water with seven-hole test section and Table II for Freon 


113 with seven-hole test section. 
Teak 
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